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Adiabatic Demagnetization Refrigerators (ADR) have many advantages that are prompting a resurgence 
in their use in spaceflight and laboratory applications. They are solid-state coolers capable of very high 
efficiency and very wide operating range. However, their low energy storage density translates to larger 
mass for a given cooling capacity than is possible with other refrigeration techniques. The interplay 
between refrigerant mass and other parameters such as magnetic field and heat transfer points in 
multi-stage ADRs gives rise to a wide parameter space for optimization. This paper first presents 
optimization strategies for single ADR stages, focusing primarily on obtaining the largest cooling capacity 
per stage mass, then discusses the optimization of multi-stage and continuous ADRs in the context of the 
coordinated heat transfer that must occur between stages. The goal for the latter is usually to obtain 
the largest cooling power per mass or volume, but there can also be many secondary objectives, such 
as limiting instantaneous heat rejection rates and producing intermediate temperatures for cooling of 
other instrument components. 
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1. Introduction 

Compared to other low-temperature refrigeration technologies, 
ADRs tend to have the highest thermodynamic efficiency, in large 
part due to the reversible nature of the magnetocaloric effect [ 1 ] 
used in the cycle. However, their need for large magnetic fields 
and the relatively low density of paramagnetic refrigerants also 
means they tend to have the lowest cooling capacity per unit mass. 
Consequently, ADR design must focus both on minimizing heat 
loads to the cold stage, and ensuring that the stage achieves the 
maximum possible cooling capacity for its mass. Complicating 
the process is that fact that a wide variety of paramagnetic refrig- 
erants can be used, and that the cooling capacity depends not only 
on the amount of refrigerant, but also on the magnetic field used to 
drive the cycle. The interdependence of these parameters means 
that many design solutions are possible that achieve the same 
cooling capacity, but with potentially large variations in total stage 
mass. In general, though, the choice of refrigerant is guided by 
relatively straightforward selection rules [2,3], and, as will be seen, 
it is possible to develop analytic expressions for the optimal 
magnetic field as a function of the ADR’s operating temperature 
range, which can then be used to compute the refrigerant mass 
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needed to obtain the desired cooling capacity. When applied to 
single ADR stages, these yield designs for which the cooling 
capacity per unit mass is maximized. 

However, this does not necessarily achieve an optimal design 
from a system perspective. For ADRs operating over a limited 
temperature range (ratio of high end and low end temperature less 
than about 20), a single stage is likely the lowest mass solution. 
When larger operating range or notably higher cooling capacities 
are needed, system mass likely can be lowered by using multi- 
stage arrangements [4] including ADRs that cool continuously 
[5], The reduced operating range of each stage and the correspond- 
ing reduction in conducted heat loads, plus the ability to use refrig- 
erants tailored to each stage’s operating range, can result in 
significant mass savings and higher overall efficiency. The latter 
can also translate to lower peak heat rejection rates, which can 
be a design driver in systems that use cryocoolers as heat sinks. 

Multi-staging, though, introduces additional operational 
parameters, including heat transfer points, hold temperatures, 
and recycling sequences. Collectively these determine the operat- 
ing range and heat loads for each stage, and, therefore, stage mass. 
The number of free parameters can be daunting, but in practice it is 
possible to develop “rules of thumb” for the parameters that gov- 
ern the interactions between stages, and separately to develop 
expressions for, for example, the optimal magnetic field of a stage 
once the operating range is established. System designs can then 
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be built up in a systematic way from a much smaller number of 
fixed parameters. 

This paper presents those strategies and guidelines for optimiz- 
ing ADR systems, starting first with individual ADR stages, and pro- 
gressing to two-stage and then multi-stage (including continuous) 
ADRs. For single stages, where the operating range and required 
cooling capacity are fixed, optimization mainly involves determin- 
ing the correct balance between magnetic field and refrigerant 
mass. For multi-stage systems, it also involves determining the 
operating range for each stage that minimizes the parasitic entropy 
generation within the ADR. As such, the optimization process not 
only identifies the lowest mass configuration, but also maximizes 
system efficiency, which has the further benefit of potentially 
relaxing requirements for the cryogenic system acting as the ADR’s 
heat sink. 

2. Optimization of a single ADR stage 

An ADR stage consists of relatively few components: a super- 
conducting magnet (producing a magnetic field, B, at maximum 
current, I max ), a salt pill containing the magnetic refrigerant, a heat 
switch to thermally link the salt pill to the heat sink, and a suspen- 
sion system that mechanically supports the salt pill within the 
magnet bore. For sensitive instruments, there may also be a need 
for magnetic shielding, either passive (ferromagnetic or supercon- 
ducting) or active (bucking coils). The stage is shown schematically 
in Fig. 1. 

A single-stage ADR is typically designed to operate between a 
fixed heat sink temperature, T sink , and a lower operating point (or 
operating range), T (ow . Fig. 2 is a plot of the entropy versus temper- 
ature for a stage using chrome potassium alum (CPA) salt, operat- 
ing at 50 mK using a 1.3 K heat sink, and a 2 T magnetic field. The 
figure assumes that the ADR stage has already been cycled to low 
temperature; if not, the cycle will start with an isothermal magne- 
tization from the point on the zero-field entropy curve that inter- 
sects the line B-C. 

At point A, the stage is cold but has no remaining cooling capac- 
ity, and must be recycled. The salt pill is thermally isolated, and is 
warmed (A-B) by adiabatic magnetization until its temperature 
reaches that of the heat sink. The heat switch is then closed, and 
magnetization continues isothermally (B-C). The entropy of the 
salt decreases as heat is transferred to the sink until full field is 
reached. The heat switch is opened, and the stage is adiabatically 
demagnetized (C-D) until the desired operating point is reached. 
The field is then slowly decreased to maintain temperature as heat 
is absorbed. When the field reaches zero, the cooling capacity is 
exhausted, and the stage is back at its starting point. 

In Fig. 2, the entropy difference between points D and A (AS) is 
the entropy capacity of the stage. The curves show specific entropy 


Heat Sink 



(Detectors, etc.) 

Fig. 1. Schematic of a single-stage ADR. 



T (K) 

Fig. 2. Ideal cycle for a single ADR stage using CPA salt, operating at 50 mK using a 
1.3K heat sink, and a 2 T magnetic field. 

(in units of the gas constant, R ) which is multiplied by the refriger- 
ant mass to get AS. The corresponding cooling capacity at low tem- 
perature is T (ow AS, and the heat rejected each cycle is T sink AS. 

The target cooling capacity is based on the heat loads that must 
be absorbed at low temperature, and the duration (or hold time). 
The heat loads must include both internal (or parasitic) sources, 
such as the heat switch, salt pill suspension and sensor wiring, 
and external sources: detectors, wiring, amplifiers, etc. It is not 
uncommon for the parasitic loads to represent the majority of 
the heat budget, and it should be understood that they tend to 
scale with the size of the salt pill, and so the design process is usu- 
ally iterative. The heat switch, in particular, has to be sized to reject 
heat to the sink within what is usually a prescribed period of time, 
and the open-state conduction will typically be a function of the 
size needed to achieve a suitable closed-state conductance. 

2 A. Selection of refrigerant 

Paramagnetic materials used as refrigerants are characterized 
by a spin quantum number J, a magnetic ion density, and an order- 
ing temperature. Magnetic ordering occurs below the temperature 
where the energy of nearest neighbor interactions is comparable to 
the thermal energy. The interaction competes with the effect of an 
applied magnetic field, reducing the latter’s ability to affect the 
entropy and produce cooling. Effectively the ordering temperature 
is a lower limit to the useful range of a paramagnetic refrigerant. 

The ordering temperature is positively correlated with J and the 
magnetic ion density. Consequently, the choice of refrigerant for an 
ADR stage is driven by the lowest temperature at which the stage 
needs to operate. Lower temperature requires the use of less dense, 
lower j-value materials, with correspondingly lower refrigeration 
capacity. 

For reviews of low temperature refrigerants, their cooling capa- 
bilities and selection guidelines, see [2,3], Here we note that the goal 
is typically to obtain the highest cooling capacity per volume in the 
salt pill, although some consideration must be given to the chemical 
properties and, hence, how the material must be packaged. 

Low temperature materials (for use below, say, 100 mK) are 
almost exclusively hydrated salts [6]. In order to make thermal 
connection between the salt and other components, the salt is 
grown from solution inside a hermetic container, onto a thermal 
bus. The bus is constructed from thermally conductive [7-10] 
material, such as gold or copper wires. Some salts, like ferric 
ammonium alum (FAA), are highly corrosive to copper, whereas 
CPA and CMN are not. This leaves gold as the only practical choice 
for FAA. Depending on the extent of the bus needed, its cost or 


142 


P. Shirron / Cryogenics 62 (2014) 140-149 


labor-intensive fabrication time may be factors which bias the 
design away from a thermodynamically superior material to one 
that can be more easily (or inexpensively) utilized, in any event, 
the volume occupied by the salt pill container and thermal bus 
reduce the amount of salt that can be accommodated in the mag- 
net bore, and this inefficiency must be weighed against the benefit 
of materials (if any exist) that do not require such structures. 

Higher temperature materials, such as gadolinium gallium gar- 
net (GGG), tend to be available as single-crystals, and stable in any 
environment. Most, in the temperature range of use, have high 
enough thermal conductance and low enough boundary resistance 
to not require an integral thermal bus. Whether one is used or not 
depends on the specific material’s properties and the heat fluxes 
that must be supported. That is, a thermal bus may be advanta- 
geous if the thermal gradients within the salt pill reduce the real- 
izable cooling capacity by more than the reduction in salt volume 
from its inclusion would. 

2.2. Optimization of magnetic field and refrigerant mass 



B (Tesla) 

Fig. 3. Specific entropy of CPA at 1.3 K. The zero-field entropy for CPA at 50 ml< is 
indicated. 


The cooling capacity of an ADR stage operating is determined by 
both the amount of salt used and the peak magnetic field used to 
drive the cycle. Since a reduction in one can be compensated for 
by an increase in the other, and since the magnet (and shield) will 
typically dominate the ADR’s mass, there are potentially significant 
mass savings to be gained by correctly balancing the salt mass and 
magnetic field. 

The process involves developing a relationship between salt 
mass and peak field that yields the desired cooling capacity, then 
identifying the combination having the lowest mass. As shown in 
Fig. 2, cooling capacity is proportional to the difference between 
the entropy of the salt at zero field and the operating temperature, 
SCO, T/ow). and at peak field and the sink temperature, S(B, T sink ). 

Entropies can be calculated using the function: 

s(B. T)/R = xcoth(x) - (2/ + 1 )xcoth((2/ + l)x) 


x EE li B gB eff /2k B T = (0.336 K/T )gB eff /T 

(where g is the Lande g-factor, is the Bohr magneton, and k B is 
Boltzmann’s constant), using the internal field approach to model- 
ing the zero-field entropy [4]: 


Beff = 



( 2 ) 


where b is the background field generated by each ion’s nearest 
neighbors. 

Assuming the operating and sink temperatures are fixed, the 
entropy capacity is then only a function of the peak field. As an 
example, we consider the case of a CPA stage operating between 
1.3 K and 50 mK. Fig. 3 shows the specific entropy of CPA at 1.3 I< 
as a function of peak magnetic field. 

if the required cooling capacity is AQ, the required salt mass (in 
moles) as a function of magnetic field is 


tlsalt(B) 


AQ/T tow 

S(B, T sink ) - S(0, T/ ow ) 


( 3 ) 


The result is shown in Fig. 4 for the example case of CPA oper- 
ating at 50 mK under a 1 pW load for 24 h. 

Translating the salt mass into a magnet bore volume, and from 
there a magnet mass, requires some details of the salt pill design, 
such as the aspect ratio (length to diameter), the fraction of the 
volume that will be devoted to the thermal bus and hermetic struc- 
ture, and the gap between the salt pill and magnet bore. For many 



Fig. 4. Moles of CPA required to absorb 1 gW for 24 h at 50 mK, when demagne- 
tized from 1.3 I< and magnetic field, B. Ideal (Carnot) efficiency is assumed. 


of these, general guidelines will suffice. The salt typically occupies 
at least 90% of the salt pill volume, and a gap of 0.75 mm between 
the salt pill and magnet bore is adequate, given the rigidity of typ- 
ical suspension systems. 

The issue of aspect ratio is more complex, as it is a key param- 
eter affecting the thermal conductance of the salt pill’s thermal bus 
[7,10], and is one of several parameters (along with magnet wind- 
ing length, wire diameter and magnetic shield configuration) that 
affect the spatial distribution of the magnetic field in the volume 
of the salt, and hence the salt’s overall cooling capacity. 

For simple solenoids, the field decreases rapidly toward the end 
of the windings. For a semi-infinite solenoid, the field at the open 
end is half the central field. For smaller aspect ratio magnets, the 
same winding density produces a lower central field, but the field 
is more uniform along the length. Adding ferromagnetic shielding 
will both increase the central field, and further improve axial uni- 
formity. Superconducting shields and bucking coils, on the other 
hand, tend to decrease the central field and increase gradients. 
Their advantage is that the shield components are less massive 
than ferromagnetic shielding. But as the magnet size must be 
increased to compensate for the lower average field, there is gen- 
erally no net mass advantage. 

As a result, it has become common to use ferromagnetic shields 
to provide basic control of fringing fields in ADRs. Some basic 
design rules were developed by Hagmann and Richards [14]. The 
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shield is sized so that when the flux returned outside of the wind- 
ings in concentrated into the cross section of the shield, the shield 
just reaches saturation. The flux can be calculated from the central 
field, and the inner and outer diameters of the windings. This 
would seem to give a simple means of estimating shield mass for 
an ADR stage, and that is true if the shield is a simple cylindrical 
shield. However, the mass of a shield can be reduced significantly 
by contouring the shape to reflect the field profile along the axis 
and at the ends. An optimization of this sort requires specialized 
magnetic field modeling software, such as MAGNETO [15], to com- 
pute fields for asymmetric geometries. 

Consequently, it is beyond the scope of a general optimization 
analysis to include magnetic shields. We note, though, that for a 
given central field, ferromagnetic shield thickness - and mass - 
is minimized when the magnet’s outer diameter is minimized. 
Thus when optimizing aspect ratio and magnetic field, some 
weight could be given to the combination of parameters that yields 
the smallest magnet diameter. In fact, when the optimal aspect 
ratio and magnetic field are based on minimum stage mass, the 
resulting magnet diameter is the minimum. 

2.3. Optimization of aspect ratio 

To demonstrate the effect of varying aspect ratio, we choose a 
point design that is representative of low temperature ADRs in 
terms of operating temperature and cooling capacity, in this case 
a salt pill containing 0.5 mol of CPA, and magnet with an average 
field of 2 T. For a given aspect ratio, one can compute the salt pill’s 
length and diameter, and the magnet’s length and inner diameter. 
(It is assumed that the two lengths are equal.) From these, the 
outer diameter of the windings needed to achieve the required 
field (2 T) is computed. Mathematica™ [11] is used to perform 
numerical integration of Maxwell’s equations for a loop current 
over the extent of the magnet’s windings and the salt pill volume. 
Magnet mass is computed using an average density of 5 g/cm 3 . 

One final input to the calculation is the engineering current 
density, j E , of the superconducting wire, which itself is dependent 
on the maximum field that is produced in the windings. In the pro- 
duction of ADRs for space applications, small diameter magnet 
wire is typically used in order to limit the current required. The 
peak current affects both the power draw of the control electronics, 
and static conduction and ohmic dissipation in the dewar’s current 
leads. The latter are often one of the major heat loads on the cryo- 
genic system. Optimization of the current therefore involves trades 
between the spacecraft power and thermal systems, the control 
electronics, the instrument, and the dewar, as well as such consid- 
erations as the commercial availability of small diameter wire. 

Magnets for space ADRs have tended to use currents in the 
range of 2-5 A. The ADRs built for XRS 12] on Astro-E, -E2 and 
for SXS on Astro-H [13] all used 0.1 mm diameter multifilamentary 
NbTi wire, and a maximum current of 2 A to produce fields in the 
range of 2-3 T. The corresponding j E is about 2.0 x 10~ 4 A/cm 2 . 
Representative critical current values (at 4.2 K) for 15S40 NbTi 
wire supplied by Supercon, Inc. [16], used in the Astro-H/SXS mag- 
nets, is shown in Fig. 5. 

Note that at 2 A, the wire is superconducting at fields greater 
than 7 T. While laboratory magnets are routinely operated to 
within a few % of their critical current, for space instruments, some 
margin is usually added (the statistics here are very limited). For 
SXS on Astro-H, the required margin on current and field was 
25%. In reality both the 2 T and 3 T magnets could sustain currents 
of in excess of 2.7 A before quenching. 

Magnet mass versus aspect ratio is shown in Fig. 6. While the 
mass has a minimum at about 2.75, the minimum is rather broad 
and shallow, and suggests that the aspect ratio can be adjusted 
as needed for other concerns (such as meeting diameter or height 



Fig. 5. Critical current versus magnetic field for multifilamentary NbTi wire (15S40) 
supplied by Supercon, Inc. 



Aspect ratio 

Fig. 6. Mass of a 2 T (average field) magnet as a function of salt pill aspect ratio for a 
0.5 mol CPA salt pill. 


constraints, or optimizing the thermal bus) without significant 
penalty in performance. 

2.4. Optimization of magnetic field 

Using 2.75 for the salt pill’s aspect ratio, Fig. 7a presents the 
estimated mass of both the salt pill and magnet for an ADR stage 
using CPA that is sized for a cooling capacity of 1 pW for 24 h at 
50 mK, and assuming the salt pill is demagnetized from T high = 1.3 - 
K. Fig. 7b shows the same analysis for GGG sized for 20 pW for 24 h 
at 1 K using a 4.25 K heat sink. The heat loads chosen represent the 
same entropy absorption for both stages, but in practice the cool- 
ing capacity does not appreciably affect the optimal field value. 

At low fields, a very large salt mass is needed due to the rapidly 
diminishing entropy capacity. At high fields, the salt’s entropy is 
already suppressed close to zero, and additional field does not 
appreciably increase the cooling capacity. 

A fuller mapping of the optimal magnetic field provides a useful 
tool for optimizing stages for any combination of operating param- 
eters. Fig. 8 shows the results for CPA over the range of T (ow = 50- 
300 mK, and T high = 0.3-2. 6 K. It turns out the two dependences 
are separable, as the function: 

Bo Pt imai CPA = 0.572 TZ 35 ^ (in Tesla) 


( 4 ) 
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Fig. 7. Component masses for (a) a CPA stage designed to operate at 50 mK from a 1.3 K heat sink, and (b) a GGG stage designed to operate at 1 K from a 4.2 K heat sink, as a 
function of magnetic field. 




Fig. 8. Optimal magnetic field for a CPA stage as a function of demagnetization temperature and operating temperature. 


The largest deviations are at the high end of the range for sink 
temperature, where the optimal field is overestimated by less than 
3%. 

A similar analysis for FAA yields: 

Bo P uma, FM = 0.355 C 45 C 2 (in Tesla) (5) 

As an example, consider the ADR flown with the XRS instru- 
ment on Astro-E and Astro-E2, for which Thigh ~ 1.5 K (using a 
1.3 K helium bath as a sink), and T low = 0.06K. From Eq. (5), the 
optimal field is 1.69 T. In practice, since the optimal design point 
is at a very shallow minimum, somewhat larger magnetic fields 
(2 T for XRS) are often chosen to provide margin on operating con- 
ditions on orbit. 

3. Optimization of two-stage ADRs 

Two-stage ADRs that are configured in parallel are essentially 
independent stages that use the same heat sink, and are optimized 
separately as for single stages. When the stages are arranged in ser- 
ies, as in Fig. 9, the upper stage (S2) has two roles: to absorb heat 
from the lower stage (SI ) in order to regenerate its entropy capac- 
ity, and to absorb heat loads during the hold time (and thereby 
reduce the heat load on SI). The transfer point between stages, 
Tj, and S2’s hold temperature, T 2 , are parameters that need to be 
optimized. 

Recycling involves two separate heat transfers, from SI to S2, 
and from S2 to the heat sink, not necessarily in that order. In fact, 
there are many possible sequences, including 


Heat sink, T sink 



Fig. 9. 2-stage ADR coupled in series. The two stages are stationed at T 2 and T (rm . 
during the hold time, and heat is transferred between stages at 7,. 


1. Recycle SI; recycle S2; demagnetize SI and S2. 

2. Recycle S2; recycle SI; demagnetize SI and S2. 

3. Recycle S2; recycle SI; recycle S2; demagnetize SI and S2. 
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Schematically these sequences are shown in Fig. 10. Each recy- 
cling sequence begins at the end of a hold time, and each recycling 
operation follows that of a single stage. For simplicity, the heat 
transfer point between SI and S2 (T,) is the same as S2’s hold tem- 
perature ( T 2 ). In all sequences, SI and S2 absorb heat during the 
hold time. In sequence 1, S2 then absorbs entropy from SI, and 
finally rejects all entropy absorbed during the cycle to the sink. 
In sequence 2, the order of heat rejection is simply reversed. In 
sequence 3, the first recycling of S2 rejects the entropy it absorbed 
during the hold time, and the second recycling rejects the entropy 
absorbed from SI. By recycling S2 multiple times, the stage can be 
smaller, sized for the larger of the two entropy loads, not the sum. 
The trade-off is a somewhat longer recycle time. Note, though, that 
each recycle of S2 is shorter than for sequences 1 and 2 since less 
heat is transferred each time. 

The choice between sequences 1 and 2 essentially depends on 
whether T, is greater than or less than T 2 . The largest entropy 
change for the cycle will be obtained when S2 completes its cycle 
at the higher of the two temperatures. To illustrate the point, 
Fig. 1 1 shows the entropy diagram for the upper stage of a 2-stage 
ADR designed for a 1.8 I< heat sink, using T, = 0.8 K and T 2 = 0.5 K. 
Using sequence 1 , after S2 is recycled, it is cooled to 0.5 K for the 
hold period. At the end of the hold, the temperature is raised to 
0.8 K where heat is absorbed from SI, then S2 is recycled again. 
The entropy capacity of the stage is much larger than if the 
sequence were reversed. 

3.3. Optimization o/T, and T 2 

Depending on how the ADR and the components it cools are 
configured, the heat loads on SI will be dependent either on the 
heat sink temperature, or on the intermediate temperature, T 2 . 
The load (detectors, amplifiers, etc.) and salt pill suspension gener- 
ally have direct links to the heat sink, but some designs will include 
a heat intercept to S2. In the latter case, S2 will absorb the majority 
of the heat conducted through those components, and the smaller 
heat load reaching SI will be dependent on T 2 . The conducted load 
from heat switch HS1 will necessarily depend on T 2 . 

Entropy absorbed by SI during the hold time is transferred to 
S2, before rejecting it to the sink. The heat load on SI therefore 
affects the required cooling capacity of both stages. The required 
entropy capacity of SI is 

AS, = QwM + Qs2(T2) AWCi (6) 

' low 

where Qi sink (T S i„k) is the conduction to SI through components 
linked to the heat sink, Qs 2 (T 2 ) is the conduction through 


components linked to S2, At hold is the hold time, and c, is the heat 
absorption efficiency (the “usable” cooling capacity as a fraction 
of ideal). For stages operating below 200 mK, the efficiency tends 
to be in the range of 70-80%, where warmer stages can achieve effi- 
ciencies of 80-90%. 

The required entropy capacity for S2 is then 

^At M )/s 2 (7) 

where Ch^ff's-ink) is the conduction to S2 through components 
linked to the heat sink, e 2 is the efficiency of stage 2, and £,_> 2 is 
the efficiency of the entropy transfer from SI to S2. (1 - e,^ 2 ) is 
essentially the ratio of the temperature gradient across HS1 divided 
by the average temperature of SI and S2 as heat is transferred. Typ- 
ically this will be 90%, although larger or smaller gradients can be 
used depending on the heat switch conductance. 

The above expression assumes that either recycling sequence 1 
or 2 is being used, for which the capacity of S2 is based on entropy 
absorption from SI and during the hold time. When sequence 3 is 
used, the expression should include only the larger of the two 
entropy terms. 

The entropy capacity of both stages can also be expressed in 
terms of their cycle parameters: 

AS, = n sa(t| (s(0, T low ) - s(B, , T,)) (8) 

AS 2 = n sa it 2 (s(0, T,) - s(B 2 ,Tsink)) (9) 

where n sa i h and n sa i t2 are the salt masses used in SI and S2, and B, 
and B 2 are the peak fields. In the expression for AS 2 it is assumed 
that T, is larger than T 2 . If not, the zero-field entropy term should 
be replaced with s(0, T 2 ). 

The effects of varying T, and T 2 can now be seen. Decreasing T 2 
lowers the heat load on SI, and the subsequent heat transfer to S2, 
but it increases the entropy generation associated with the heat 
load on S2 during the hold time. Below the optimal value of T 2 , 
the increase in salt mass and magnetic field for S2 needed to 
accommodate the elevated entropy generation rate more than off- 
sets the reduction in salt mass and magnetic field for SI and S2 
associated with the lower heat load on SI. 

Decreasing T, increases the entropy capacity of SI, but decreases 
the entropy capacity of S2. The optimal value of T,- is the point 
where the combined mass of the salt pills and magnets for SI 
and S2 is a minimum. 


as 2 = (as,/£,^ 2 + 


Sequence 1 


Sequence 2 


Sequence 3 



Time (arb. units) 


Time (arb. units) 


Time (arb. units) 


Fig. 10. Typical recycling sequences for a 2-stage ADR. Recycling begins at the end of the hold time (time equal zero). 
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Fig. 11. Entropy diagram for the upper stage of a 2-stage series ADR, depicting sequence 1 (left) and sequence 2 (right) for the case when T, is greater than T 2 . The diagram is 
based on the properties of gadolinium lithium fluoride (GLF) at representative fields. 


4. Optimization of multi-stage and continuous ADRs 

In the preceding sections, it has been assumed that the number 
of ADR stages is fixed, and the parameters affecting their operation 
are varied in order to obtain the highest cooling capacity per unit 
mass. But for an ADR using stages connected in series between 
two temperatures, the number of stages can also be considered 
an adjustable parameter. Multi-staging can reduce the overall mass 
by (1 ) reducing the magnetic field requirements for each stage, and 
(2) reducing parasitic heat loads and associated entropy generation. 

A simple case to examine is the use of a one or two ADR stages 
to span an increasing temperature range. Using 50 ml< as a “stan- 
dard” operating temperature, and a fixed cooling capacity of 
0.0864 J (1 pW absorbed over 24 h), Fig. 12 compares the mass of 
single- and two-stage ADRs as a function of T sink . CPA is used for 
all stages. For the two-stage ADR, the hold temperature, T 2 , and 
the intermediate temperature, T„ are both equal to the logarithmic 
average of 50 ml< and T sink . (In this case, one can use either recy- 
cling sequence 1 or 2, as they are equivalent.) The optimal mag- 
netic fields are computed from Eq. (4). 

For fields above ~2 T, the single-stage mass has a temperature 
dependence of ~ T^ k . The two-stage ADR uses much lower mag- 
netic fields in each stage, with correspondingly lower magnet 
masses. As a result, it is advantageous to transition to two stages 
when the ratio of temperatures ( T sink ITi ow ) exceeds ~30, at which 
point each stage has an operating range of about a factor of ~5.5 
in temperature. 



The advantage of multi-staging is even more pronounced when 
the effect on cooling capacity is considered. Components, such as 
heat switches, that link to an intermediate stage rather than 
directly to the heat sink will impose lower heat loads on the lower 
stage by a factor of 

c. rh / rT sink / T. \ m +' 

Y= k(T)dT / / k(T)dT* (10) 

Jt, ow / JT bw VsinkJ 

where thermal conductance is taken to be of the form k~T m , and 
T™ w -c T™, T™ nk . If TilT sink — 1/5 and m = 1, the reduction is a factor 
of 25. 

The total heat load will, in general, be reduced by a smaller frac- 
tion due to contributions that are not affected by the intermediate 
stage, and a full optimization needs to be performed for each stage 
as outlined in the previous section. However, the combined bene- 
fits of smaller magnetic fields and lower heat loads suggest that an 
ADR’s performance can be improved if, by adding a stage, the tem- 
perature span {T h i gh ITi ow ) of each stage is reduced into the range of 
5-6. 

4.1. Continuous ADRs 

The physical arrangement of stages in a series CADR [5] is the 
same as for single-shot, multi-stage ADRs, and each stage operates 
on the same cycle: giving off heat during an isothermal magnetiza- 
tion phase at “high" temperature with its heat switch closed, and 



Fig. 12. Comparison (left) of the masses of a single-stage and a two-stage ADR spanning 50 ml< to Tsink, and the optimal magnetic fields (right) for both ADRs. For the 2-stage 
ADR, T 2 = Tj = ^7(50 mK)T sink , and recycling sequence 1 is used. 
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absorbing heat at “low" temperature during an isothermal demag- 
netization phase with its heat switch open. If there is a lower stage, 
the heat absorption is done with that stage’s heat switch closed in 
order to regenerate its cooling capacity. 

The first notable difference between CADRs and single-shot 
ADRs involves the control of the lowest stage in the chain (the 
“continuous” stage). Its low and high temperatures are the same, 
and it alternates between isothermal magnetization and isother- 
mal demagnetization as needed to maintain a constant tempera- 
ture. The second difference is the duty cycle of each stage. In a 
single-shot ADR, the cascade of heat to the sink is performed over 
a small fraction of the total cycle, so that the majority of time is 
spent in the isothermal demagnetization phase. In a CADR, the 
cycling is continually repeated, with each stage spending about 
half the time rejecting heat, and half the time absorbing heat. 

This has an impact on the requirements for the heat switches, 
and consequently, for the operating range of each stage. Heat 
switches are designed to have high enough thermal conductance 
in the closed state to recycle the ADR within an allocated time, 
and low thermal conduction in the open state to minimize parasitic 
heat loads. The two requirements tend to drive the design in oppo- 
site directions - larger structures needed to obtain higher “on” 
conductance generally have higher heat leak in the “off’ state. 

The linkage is not absolute, of course. For superconducting heat 
switches, the A/L of the superconductor dictates both the on and 
off state conductances [17], But for gas-gap heat switches, the on 
conductance is proportional to the volume occupied by the con- 
ductive elements, while the off conduction is a function of the 
aspect ratio of that volume. Thus depending on the type of switch 
used, there may be design choices that allow higher on conduc- 
tance without sacrificing off conduction. 

The requirements for heat switch performance can be derived 
for different ADR configurations. Following an analysis by Kittel 
[IS, we evaluate the entropy flow of a heat switch in both the 
on and off states. During the isothermal demagnetization phase 
(absorbing heat at T fow while the warm end of the heat switch is 
T high ), the entropy load from the heat switch is 

S 1 ,„ s = (^) 7 L [ T ™k 0FF (T)dT (11) 

V T ' 1 low 4 t Iow 

and the entropy flow out during the isothermal magnetization 
(rejecting heat at T re j e ct) is 

S„ut =h- T f L ) y*— [ T "’ ea+AT 1 <on(T) dT (12) 

V ~ ' 1 re i ea JT„ jea 

where is the fraction of time that the stage is cold (also called 
the duty cycle), AT is the gradient across the switch when heat is 
being rejected, and k 0N (T) and IcoffCO are the on and off state con- 
ductances. At low temperature, the latter can be expressed as fc 0 iv(- 
T) = aT 1 and k 0FF (T) = fiT. 

In a practical system, the stage is designed to lift not only its 
own parasitic heat, but also heat from a load, which adds an 
entropy term Q/T, ow to S,-„. Assuming T™ w « TZ gh and AT « T re jec[ , 
the entropies reduce to 


Sin 


Q. f7-hold\ 1 /I jm+1 

how W /T^m+T 


(13) 


S 


out — 



jn-\ at 
J - 1 reject LA1 


(14) 


an efficiency factor that accounts for gradients within the salt pill, 
eddy current losses, etc. 

Applying these expressions to the continuous stage: T re j ect = 
T( 0W , and t/, oW /t ~ 0.5. The on and off state conductivities will 
depend on the particular switch technology used. For a stage that 
will operate below about 200 mK, the choice is limited to those 
that use metallic conductors: superconducting, magnetoresistive, 
or mechanical. For this discussion, a superconducting switch is 
assumed. The off thermal conductance is governed by the popula- 
tion of quasiparticles, giving a relatively strong temperature 
dependence. Lead (Pb) switches developed at NASA/GSFC [17] 
have /coff{T) = 2.15 T 3 ' 21 (mW/K), In the on state, k O N(T) = 110 T 
(mW/K). 

Inserting these values gives the requirement 
S„ ut = 55 AT(mW/K) > S in /s 

= ^(q + 0.256 T£ 2 ’(mW/I<)) (15) 

For a CADR targeting Q = 5 pW of cooling power at T (ow = 50 mK 
(using AT = 5 mK), and using s = 0.7, T high is constrained to be 
<0.38 K. At this temperature, the time average heat load on the 
continuous stage is 9.36 pW. 

If the size of the superconducting element is changed, for exam- 
ple to a larger A/L, both the on and off conductances change by the 
same proportion. The effect will be to make the fixed heat inputs 
(Q) smaller relative to the heat rejection capability. On the surface 
this would provide leeway to increase T Mgh , but higher heat flows 
can lead to less efficient transfers, and such adjustments should 
be approached cautiously. 

For the upper stages of the CADR, switching ratios are less of a 
driver for operating temperature, and there is greater flexibility in 
the choice of technologies. Gas-gap switches can be used at tem- 
peratures above about 0.3 K, with very high intrinsic on/off ratios. 
On state conductances on the order of k ON (T) ~ 25T (mW/K) are 
easily achieved, with off conductances on the order of /<off{- 
T) ~ 0.5 T 2 ' 4 (pW/K) 

Using these parameters, along with T/, oM /t ~ 0.5, Eqs. (13) and 
(14) reduce to 

O T 34 

Sin =r^+ 0.074^ (pW/K) (16) 

J low * low 

S out = 12.5AT(mW/K) (17) 

The direct heat load term is retained, and represents the 
entropy that each stage must absorb from the next lower stage. 
For the second stage of the CADR (connected to the continuous 
stage via the superconducting heat switch), this term is 9.38 pW 
divided by the efficiency of 0.7. The remaining parameters are: 
T/ ow = 45 mK, T high = 0.38 K, and AT = 0.038 K. The two entropy 
terms contributing to S in are 0.30 mW/K and 0.030 mW/K, the 
smaller attributable to the gas-gap heat switch. S out , in this case, 
is 0.48 mW/K - more than adequate for cascading heat to the next 
(3rd) stage. 

Since heat switch performance is not the driver for the operat- 
ing range of the 3rd and higher stages, their optimization will be 
driven by the individual and overall stage masses as described at 
the beginning of this section. The main result is that each stage 
is optimally designed to operate over about a factor of 5 in temper- 
ature. Table 1 summarizes the operating parameters for a 4-stage 
CADR [5] designed to operate between 50 mK and 5 K. 


In an idealized ADR, the heat switch must have characteristics 
that yield S out > S in . In a practical system, the realizable cooling 
capacity is typically 70-90% of ideal, depending on temperature. 
A more realistic requirement is therefore S out > S in /e, where e is 


4.2. CADR cooling power 

A final consideration for CADR design is the relationship 
between system parameters (heat switch characteristics, operating 
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Table 1 

Operating Parameters for a 4-stage CADR operating between 50 mK and 5 IC 


Stage 

Refrigerant 

B( T) 

Thigh (K) 

Tlow (K) 

Mass (kg) 

1 

60 g CPA 

0.1 


0.05 

0.5 

2 

100 g CPA 

0.5 

0.3 

0.045 

1.7 

3 

100 g CPA 

1.5 

1.3 

0.25 

2.4 

4 

65 g GLF 

4.0 

5 

1.15 

3.1 


temperatures, salt pill size, etc.) and cooling power. Unlike single- 
shot ADRs where salt pill size is dictated by heat loads and required 
hold times, CADRs do not have a hold time requirement. Instead, 
the relevant time frame for computing cooling capacity is the cycle 
period. The cycle period, however, depends largely on the time 
required to transfer the cooling capacity of the continuous stage 
to the 2nd stage. 

As a result, the size of a CADR's salt pills have far less impact on 
cooling power than the conductance of the heat switches. And in 
this regard, the continuous stage’s heat switch is most important, 
since it has a direct effect on the entropy load that must be cascaded 
through the CADR. Its on conductance can be increased with a lar- 
ger superconducting element, but this also raises the off conduc- 
tance. An effectively larger switching ratio can be obtained by 
reducing the peak temperature of the 2nd stage. In the 4-stage 
CADR, T h i gh for the superconducting switch was limited to 0.3 K, 
reducing its average load on the continuous stage from 4.2 pW to 
1 .6 pW - nearly a 30% reduction in the total heat load. 

CADR performance can also be improved by reducing the frac- 
tion of the cycle time spent in “overhead" operations, including 
ramping the upper stages over their temperature range, and open- 
ing and closing heat switches. Because the CADR magnets are 
small, the stages can easily be ramped between operating points 
within a minute (each way), and heat switch transition times can 
be limited to zero in the case of passive gas-gap heat switches, or 
to less than a minute for active versions. At these (de)magnetiza- 
tion rates, eddy currents [9] are a concern, and the salt pill’s are 
specifically designed to minimize dissipation. Primarily this con- 
sists of minimizing closed loops of high purity copper used in the 
thermal buses, and locating bulk materials used for thermal inter- 
faces outside the magnetic shields. A worst-case estimate is that 
each salt pill’s cooling capacity must be increased by 5% to account 
for eddy current losses over a cycle. 

In practice the overhead operations add 2-4 min to the CADR 
cycle. For the system described in [19], the resulting cycle time 
was ~20 min. Since the overhead operations are relatively fixed 
in duration, some improvement in performance can be realized 
by increasing the cooling capacity of each salt pill, and increasing 
the duty cycle for heat transfer. 


5. Summary 

ADRs are unique among ultra-low temperature refrigerators in 
that they can use a wide range of (magnetic) refrigerants, and their 
cooling capacity is determined both by the amount of refrigerant 
and the magnetic field used. In addition, there are many architec- 
tures that can be employed - from single-stage units, to multi- 
stage single-shot systems, to multi-stage continuous versions. 
There is, then, a very wide parameter space for optimizing their 
design and performance that can at first be daunting. 

However, a given set of operating requirements can be 
translated into an optimized system by design guidelines outlined 
in this and other references [3, 7, 9, A, 17]. These optimization rules 
developed in this paper, based on obtaining the lowest stage mass 
for a required cooling capacity, can be summarized as follows. 


• The optimal refrigerant will have an ordering temperature 
slightly below the lowest desired operating temperature. CPA 
and FAA are typical for stages operating down to 30 mK and 
50 ml<; above 0.5 K, GLF and GGG become usable. 

• The use of ferromagnetic shielding enhances a magnet’s central 
field and reduces axial and radial field gradients; superconduc- 
ting shielding and bucking coils have the opposite effect, and 
require scaling up the central field. 

• The optimal salt pill aspect ratio is ~2.75, although the mass 
minimum for this parameter is broad and shallow. Different 
aspect ratios can be used to accommodate size constraints or 
to improve thermal bus conductance without significant mass 
penalty. 

• For CPA and FAA, an expression for the optimal magnetic field is 
presented in terms of the high and low operating temperatures 
of the stage. From these, a stage operating between 50 mK and 
1.5 K would optimally use a ~2 T field. 

• ADR stage mass grows rapidly as the required field exceeds 
~2 T, and/or the operating range exceeds a factor of ~30. 
Increasing the number of stages reduces mass by reducing the 
required magnetic fields and cooling capacities. The intermedi- 
ate temperature of a 2-stage ADR should be chosen to minimize 
the entropy absorption during the hold time, and will be close 
to the logarithmic average of the low end and heat sink 
temperatures. 

• For multi-stage ADRs, the recycling sequence has a significant 
affect on the mass of each stage. Optimization is based on the 
relative size of internal (parasitic) and external heat loads. Recy- 
cling upper stages multiple times can significantly reduce their 
required cooling capacity. 

• A continuous ADR (CADR) has significantly lower mass than a 
single-shot ADR for the same cooling power. Performance is dic- 
tated by cycle period, and is far more driven by heat switch con- 
ductances than salt pill size. 

• The operating range of each stage in a CADR is dictated by the 
on- and off-state conductances of the lowest temperature heat 
switches, and is optimal at r W g h /T tow ~ 5. 
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